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ABSTRACT 


Research  progress  has  been  Bade  toward  two  goals  in  the  first  contract 
year,  patient  studlsa  and  aodel  studies.  In  the  patient  studies  we 
investigated  the  relationship  between  assn  airway  pressure  and  lung  voluae 
during  high  frequency  low  tidal  voluae  ventilation  (HFT).  Patients 
requiring  neohanleal  ventilatory  support  for  treatment  of  respiratory 
insufficiency  were  studied  by  laposlng  rapid  (1-10  He)  osoillatlona  with 
low  tidal  voluwes  (50-150  al)  at  a  constant  naan  airway  pressure  of  5  oa 


HgO.  /Even  though  aean  airway  pressure  was  oonstant,  lung  voluae 
Increased  substantially  during  the  osolllation  period  in  7  of  8  subjects 
ss  Indicated  both  by  an  increase  in  tbcraoo-abdcalnal  dlaenslona  and  by  an 
Increase  in  respiratory  systaa  relaxation  pressures  after  the  oscillations 
wens  stopped.  For  each  patient  in  whoa  these  changes  occurred,  the  degree 
of  flung  inflation  rose  progressively  with  increases  in  either  frequency  or 
tidal  voluae.  These  results  suggest  that  lung  aeohanics  substantially 
influence  the  lung  voluaea  that  are  achieved,  even  at  fixed  airway 
pressures,  during  HFV.  The  iaplloations  of  these  findings  to  ventilatory 
const  casualties  are  discussed. 


vIn  the  aodel  studies  the  pressure  drop  during  sinusoidal  aean  flows  in 
a  four  generation  network  of  rigid,  unifora  dlaaster,  syuaetrlcally 
branching  tubes  was  studied.  The  data  Obtained  were  analysed  via  a 
process  of  Fourier  decoaposltlon.  The  results  showed  that  the  pressure 
slggils  consist  aalnly tof  a  doa Inant  component  at  the  excitation  frequency 
("fuMaaental'T'and  a  "first  haraonie'of  saaller  asgnitude.  He  found  the 
aagnltude  and  phase  of  the  fundaaental  to  correlate  oloaely  with  classical 
predictions  as  long  as  the  r^rsasttr  rBsynnlrtt  awake*  divided  by  Hnasrslsj 
v^janaber]^ was  less  than  200.  Afor  values  of  this  paraaeter  greater  than  200, 
the  observed  pressure  drops)  were  considerably  higher.  He  have  attributed 
this  change  in  behavior  to  the  onset  of  turbulence  in  the  branching 
network.  These  results  were  eaployed  in  a  theoretical  aodel  to  predict 
total  airway  realatanoe  of  nhe  dog  lung.  The  predictions  were  co spared  to 
physiological  asasureaents  tar  others  and  were  found  to  be  in  excellent 
agreeaent.  These  aodels  aay  be  of  value  in  predicting  the  pressure- drops 
that  will  prevail  under  HFT  circuast voces. 


Foreword 


a.  Figures  1-16  bare  been  submitted  for  publication  at  this  tlae  and 
are  not  protected  by  copyright. 

e.  For  the  protection  of  busan  subjects  the  inTestlgator(s)  have 
adhered  to  the  policies  of  applieabla  Federal  Law  45CFM6. 
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I.  UTTXODOCTION 

The  garni  of  this  oon tract  is  to  understood  the  lntoraotions  between 
pulnonary  ■school os  sad  gas  exchange  under  the  conditions  of  low  vol use- 
high  frequency  ventilation  (HFV).  In  the  first  oon tract  year  we  have  node 
progress  toward  this  coal  in  2  areas,  clinical  studies,  and  nodal 
studies.  The  Hjcr  thrust  of  the  cllnloal  studies  has  been  to  deters! ns 
tbs  relationships  anoag  applied  tidal  volume  and  frequency  end  the  seen 
lung  volune  about  which  these  oscillations  finally  occur.  Through  our  HIT 
subcontract  we  have  begun  to  develop  s  sound  fluid  dynaale  basis  for  the 
pressure-flow  relationships  that  prevail  during  HTV.  Ve  will  review  the 
progress  node  in  each  of  these  areas  separately. 


II.  CLMICAL  STUDIES:  THE  ZKPLUEXCE  OP  HFV  OR  FURCTIORAL  RESIDUAL 
CAPACITY 

A.  Introduction 

The  tidal  voluaes  used  in  the  application  of  HFV  are  snail  in 
cooper ison  to  those  used  in  conventional  neohanlcel  ventilation;  as  a 
result  It  has  been  assuaed  that  the  average  alveolar  pressuree  during  R IT 
will  be  lower  than  those  during  conventional  nechanical  ventilation. 

These  presuned  lower  pressures  underlie  nany  of  the  potential  advantages 
of  HFV.  For  exanple,  lower  alveolar  pressures  night  reduce  the  frequency 
of  pmuno thorax  and  other  types  of  barotrauaa  associated  with  nechanical 
ventilation.  Seooadly,  the  lower  pressures  nay  ninlalse  the  adverse 
cardiovascular  effects  of  nechanical  ventilation,  particularly  the  cyclic 
reduction  of  venous  return  synchronous  with  lung  inflation  during 
conventional  aechanloal  ventilation  (1). 

These  potential  advantages  of  HFV  are  based  on  the  assumption  that 
naan  airway  pressure,  alveolar  pressure  and  lung  volune  correlate  closely 
and  predictably  during  HFV.  However,  several  recent  aninal  studies  (2-h) 
have  questioned  the  validity  of  this  assumption.  Ve  exanined  the 
relationship  between  naan  airway  pressure  and  lung  volune  during  HFV  in  6 
patients  with  respiratory  insufficiency  to  provide  lnfometlon  about  the 
potential  clinical  implications  of  these  aninal  studies.  Our  results 
demonstrate  that  during  HFV  the  naan  airway  pressure  correlates  poorly 
with  lung  volune.  Ve  found  elevated  lung  volumes  during  the  application 
of  HFV  despite  aalntenanee  of  a  fixed  naan  airway  pressure. 

B.  Methods 

1 .  Patients 

Bight  patients  requiring  nechanical  ventilatory  support  served  as 
subjeots  for  this  study.  Table  1  gives  tbe  age,  sex,  weight,  diagnosis, 
static  respiratory  systen  oonplianoe  (neasured  prior  to  HFV),  and 
endotracheal  tube  else  of  aaeh  patient.  Patlaots  1-7  required  chronic 
nechanical  ventilatory  support  via  cuffed  tracheostomy  tubes.  Patient  8, 
an  asthma tio,  required  ventilatory  asaistanM  for  several  days  and  was 
ventilated  via  a  cuffed  orotracheal  tube.  This  study  was  approved  by  the 
Hunan  Subjects  Committee,  and  each  patient  or  his/her  guardian  gave 
informed  consent  prior  to  participation. 

PMVIOUS  M«t 
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2.  Experimental  Apparatus 


Tha  apparatus  tasad  Is  Ulustratad  la  Figure  1*  Tba  high  frequency 
ventilator  was  a  servo-oootrolled  linear  aagnetlc  actor  with  an  attached 
piston.  Va  used  this  to  deliver  sinusoidal  tidal  voluaea  of  50-150  nl 
over  tba  desired  frequency  range.  A  Flalsoh  #2  pnauaotaobograph,  coupled 
to  a  differential  pressure  transducer  (Yalldyne  ♦  50  oa  H2O),  calibrated 
for  each  tidal  voluae  and  frequency  used*  was  plaoed  In  the  airway  between 
the  HFT  piston  and  the  patient.  The  output  of  the  pncuaotachograph  was 
electronically  Integrated  to  deteralne  the  tidal  voluae  delivered  to  the 
patient.  The  ventilator  ciroult  included  a  weighted  10  L  water- aealed 
splroaeter  filled  with  oxygen  which  provided  gas  to  the  airway  opening  at 
a  oonstant  pressure  of  5  on  H2O.  The  splroaeter  was  oonneoted  via  a 
valve  to  the  patient's  airway.  A  rapidly  responding  transducer  (Alltech 
MS10-S)  was  used  to  aeasure  airway  pressure  at  the  junotlon  of  the 
endotracheal  tube  and  the  pnauaotaobograph  through  a  lateral  port.  The 
tubing  dlaaeter  at  this  point  was  1.7  oa.  In  patient  8,  the  splroaeter 
circuit  was  not  used  to  aalntaln  a  oonstant  airway  pressure  during  HFT. 
Instead ,  a  high  iapedance  bias  flow  as  previously  described  (5)  provided 
fresh  gas  to  the  systca  at  a  point  between  the  ventilator  and 
pneuaotaehograph.  This  systea  provided  fresh  gas  to  the  airway  and 
rcaoved  alxed  expired  gas  via  an  exit  port  in  such  a  Banner  that  aean 
airway  pressure  reaalned  fixed  at  5  oa  HpO. 

Lung  voluae  was  aonltored  with  a  respiratory  Inductance  plethysaograph 
(Sea pi trace  B)  that  was  used  to  aeasure  the  cross  sectional  area  of  the 
thorax  and  of  the  abdoaen.  This  device  was  calibrated  by  delivery  of 
known  voluaes  of  gas  to  the  respiratory  systea  of  eaoh  patient.  A 
standard  splroaeter  was  used  for  this  purpose  and  the  patients  were  in  the 
posture  assuaed  during  HFT. 

3.  Protocol 

Prior  to  the  trial  of  HFT,  each  subject  received  several  large  breaths 
of  oxygen  (1000-2000  al)  and  was  allowed  to  exhale  passively  until  a 
stable  lung  voluae  was  reached  (usually  within  10  seoonds).  The 
endotracheal  tube  was  connected  to  the  HFT  oircuit  and  the  valve  to  the 
weighted  splroaeter  was  opened.  As  a  result  of  the  5  oa  B2O  pressure, 
lung  voluae  rose  slightly  (usually  100-200  al).  After  the  lung  voluae  had 
stabilised  (2-3  see),  we  turned  on  the  oscillator  at  a  preset  frequency, 
end  adjusted  the  voluae  of  oscillation  to  deliver  the  desired  tidsl  voluae 
to  the  subjeot.  Oscillations  were  continued  until  a  stable  lung  voluae  as 
reoorded  by  the  respiratory  Inductance  plethysaograph  and  the  splroaeter 
was  reached  (usually  15-30  seconds).  Ve  then  closed  the  valve  to  the 
splroaeter,  discontinued  the  osolllstlonsi  and  aeasured  the  relaxation 
pressure  In  the  systea.  Ve  did  not  use  runs  in  whloh  a  steady  relaxation 
pressure  was  not  achieved.  This  protocol  was  repeated  to  deteralne  the 
changes  In  lung  voluae  over  a  range  of  HFT  frequencies  (1-10  Hs)  and  tidal 
voluaes  (50-150  al).  Although  arterial  blood  gas  tensions  were  not 
aeasured,  the  brief  duration  of  these  trials  (approxlaately  30  seoonds) 
would  probably  not  have  resulted  In  a  rise  In  the  tension  of  COj  of  wore 
than  2  or  3  Terr  even  If  no  gas  transport  oocurred  between  the  patient  and 
the  oxygen-filled  splroaeter.  Conventional  aeohanleal  ventilation  was 
used  to  support  the  subjects  between  trials. 


C.  Raa lilt* 


Despite  maintenance  of  a  oonstant  Man  airway  pressure,  In  7  of  6 
subjects  Man  lung  voIum  increased.  Tba  incraaaa  in  lung  voIum  a  bora 
functional  roaidual  capacity  (PRC)  aa  a  function  of  frequency  ia  given  in 
Figure  2.  Aa  the  oaoillatory  frequency  increaaed,  lung  toIum  roaa 
progressively  in  the  7  subjeots.  Further,  for  a  given  frequency  of 
oaoillatlon  tba  lung  voIum  during  oaelllation  increaaed  to  a  greater 
extent  at  higher  tidal  volunea  than  at  lower  tidal  voluMa.  The  nagnitude 
of  the  lung  voIum  incraaaa  differed  aubatantially  among  subjects. 

Subject  8  displayed  a  voIum  change  at  the  lowest  frequency  (1  Hz)  and 
tidal  voIum  (50  nl)  tested,  while  the  other  subjects  displayed  voIum 
changes  only  when  higher  frequencies  or  tidal  voIums  were  used.  Subject 
1  did  not  demonstrate  a  change  in  lung  voIum  at  any  of  the  ventilatory 
ratea  or  tidal  voIums  studied. 

Figure  3  gives  the  data  shown  in  Figure  2  aa  a  function  of  the  product 
of  frequency  and  tidal  voIum.  These  graphs  show  that  the  increase  in 
lung  voIum  during  HFT  was  related  Mrs  closely  to  the  amplitude  of 
tracheal  flow  than  to  frequency  or  tidal  voIum  alone.  Thus,  in  patients 
2,3,5,  and  7  the  change  in  ling  voIum  at  a  given  flow  was  the  sane 
regardless  of  the  frequency  and  tidal  voIum  used  to  achieve  that  flow. 

Figure  4  shows  the  relationship  between  the  respiratory  system 
relaxation  pressure  after  oscillation  and  the  lung  voIum  increment 
Masured  by  the  inductance  plethyamograph  during  each  HFV  trial.  The 
relaxation  pressure  reoorded  after  cessation  of  each  HFV  trial  correlated 
closely  with  the  change  in  lung  voIum  noted  during  the  same  trial. 

The  pressure  swings  measured  at  the  airway  opening  during  oscillation 
were  symMtrlo  around  the  Man  airway  pressure  of  5  on  HyO  during  trials 
in  which  lung  voIum  remained  constant. '  However,  in  trials  when  lung 
inflation  occurred,  the  pressure  swings  be case  asymmetric,  with  negative 
deflections  exceeding  positive  deflections  (Figure  5A).  In  Figure  5B  a 
representative  respiratory  inductance  plethysMgraphlc  recording  during 
HFV  is  shown.  This  figure  shows  the  lnstriamnt  output  during  several 
large  breaths,  followed  by  disconnection  from  the  ventilator  And 
exhalation  to  FRC.  The  lung  voIum  change  when  the  subject  was  attached 
to  the  HFV  circuit  ("a”  in  Figure  5)  is  displayed,  followed  by  a  further 
increment  in  lung  voIum  during  the  HFV  trial  Itself  (*b"  in  Figure  5). 

D.  Diaminainn 

Our  findings  dcMnstrate  that  there  can  be  a  dissociation  between  lung 
voIum  and  Man  airway  pressure  during  HFV.  Tbs  Increase  in  lung  volume 
cannot  be  explained  by  reflex  or  voluntary  respiratory  efforts,  since 
subjects  2  and  7  suffered  from  paralysing  neuroauseular  disease  that 
limited  their  vital  oapaoltles  to  less  than  100  nl,  a  voIum  far  less  than 
the  Manured  Increases  in  the  lung  voIum  during  HFV.  The  lung  volume 
ehange  cannot  be  a  measurement  artifact  slnoe  the  respiratory  system 
relaxation  pressure,  determined  after  HFV  trials  in  which  the  voIum  rose, 
exceeded  the  airway  pressure  Mintained  before  and  during  each  trial. 

Thus,  mean  airway  pressure  during  HFV  delivered  in  this  manner  is  lower 
than  mean  alveolar  pressure.  These  findings  have  important  implications 
for  the  safe  olinloal  use  of  this  approach  to  Mohanleal  ventilation. 
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A  wjor  clinical  goal  of  lev  tidal  voluw  ventilation  is  to  provide 
aoobanieal  ventilatory  support  for  patients  at  loner  nean  airway  pressures 
(6-8) ,  assuming  that  the  inoidenee  of  adveree  side  effects  will  dialnish 
if  the  airway  and  alveolar  pressures  during  wofaanloal  ventilation  can  be 
nlnlalsed.  Iapliolt  in  this  asauwptlon  is  the  belief  that  wan  airway 
pressure  sowhow  reflects  alveolar  pressure.  Our  results  suggest  that 
this  assumption  is  not  always  valid  during  HFV.  Specifically,  our 
wasurenenta  at  the  airway  opening  oonflrwd  that  both  the  weighted 
spirowter  and  the  bias  flow  system  successfully  wlntained  wan  pressure 
constant  throughout  the  trials  of  HFV.  Despite  this  constant  wan  airway 
pressure,  lung  voluw  Increased  during  HFV. 

The  system  we  used  in  thew  experiwnta  is  alnilar  to  systems  that  we 
(9)  and  others  (10,11)  have  uwd  to  deliver  HFV  experimentally  and 
clinically.  Although  we  wde  no  attwpt  to  aonitor  pulmonary  gas  exchange 
during  thew  short  runs  of  HFV,  we  used  frequency  and  tidal  voluw 
combinations  which  were  within  the  range  required  to  achieve  adequate 
ventilatory  support  in  adult  human  subjects  (9,11).  Thus,  it  seems  quite 
likely  that  there  will  be  a  dissociation  between  wan  lung  voluw  and  wan 
airway  pressure  during  clinical  applicaton  of  HFV. 

Dissociation  of  wan  airway  pressure  and  lung  voluw  during  HFV 
probably  underlies  certain  obwrvatlona  wde  by  other  researchers  in 
animal  experiments.  For  example,  Simon  and  oo- workers  found  that  airway 
relaxation  pressures  after  HFV  in  dogs  exceeded  the  wan  airway  pressure 
measured  during  HFV  (2).  Further,  Robertson  and  oo-workers  noted 
increases  in  the  functional  residual  capeolty  of  paralysed  dogs  during  HFV 
(3),  although  their  protocol  does  not  clearly  state  whether  wan  airway 
pressures  remaiwd  constant.  This  study  extends  these  obwrvatlons  to 
humans,  and  documents  that  the  extent  of  this  hyperinflation  wy  be  great 
enough  to  cause  adverse  effects,  though  none  were  noted  during  the  brief 
duration  of  thew  trials. 

Even  though  this  study  was  not  dwlgned  to  determine  the  precise 
wcbanlsns  accounting  for  dynamio  hyperinflation  during  HFV,  we  can 
suggest  potential  explanations  for  our  obwrvatlons.  Dynamio  hyper¬ 
inflation  will  occur  when  the  voluw  of  gas  entering  the  lung  during 
inspiration  exceeds  the  voluw  of  gw  leaving  the  lung  during  expiration. 
For  this  to  oeour,  the  resistance  to  expiratory  flow  must  be  high  enough 
that  the  forces  promoting  exhalation  are  inadequate  to  allow  full 
exhalation  before  the  next  breath  is  delivered.  In  this  context,  it  does 
not  wtter  whether  the  expiratory  resistance  la  due  primarily  to  airway 
narrowing  (fixed  or  dynamic)  or  to  the  resistance  of  the  rwplratory 
pathway  of  the  HFV  circuit,  or  the  patient's  own  airway.  However,  since 
our  slrway  pressure  recordings  (Figure  5A)  showed  a  wrked  asymwtry,  it 
is  likely  that  dynamio  airway  narrowing  played  a  wjor  role  in  lung 
inflation  during  high  frequency  oscillation.  Regardless  of  the  precise 
reason  that  lung  voluw  inoreaws  dynamically,  when  this  does  occur, 
elaatlo  reooll  forces  will  increase  and  airway  resistance  will  fall. 

Thus,  lung  voluw  will  in  or  saw  until  a  olrcwstanoe  is  reached  where  the 
voluw  expired  in  the  brief  tlw  allotted  equals  the  Inspired  voluw. 

This  wy  explain  why  patient  1,  who  bad  very  noaooapllaat  lungs,  did  not 
demonstrate  an  increase  in  lung  voluw;  her  reooll  pressures  were  adequate 
to  balance  expiration  and  inspiration  war  FRC  at  all  frequenolas 
applied.  Thew  airway  and  elastic  reooll  oonsideratlons  will  apply  to  all 
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types  of  aechanlcal  ventilation,  regardless  of  whether  exhalation  is  an 
active  or  a  passive  process;  thus  a  dynamic  increase  in  lung  volume  is  a 
potential  risk  of  all  types  of  rapid  ventilation,  Including  both  high 
frequency  oscillatory  ventilation  and  high  frequency  jet  ventilation. 

Our  findings  show  that  HFV  may  not  be  a  suitable  method  of  mechanical 
ventilatory  support  for  some  patients.  For  example,  three  patients 
manifested  potentially  dangerous  increments  in  lung  volume  (over  1L) 
during  these  HFV  trials.  On  the  other  hand,  some  patients  showed  little 
or  no  evidence  of  dynamic  lung  inflation,  at  least  at  modest  tracheal 
flows.  In  previous  studies  of  HFV  in  human  subjects  (9,12,13),  we  and 
others  have  demonstrated  that  eucapnic  gas  exchange  can  be  achieved  at 
tracheal  flows  in  the  range  of  5  to  15  L/min.  In  the  current  study,  5  of 
the  8  patients  had  dynamic  shifts  in  FRC  of  less  than  500  ml  at  these 
flows.  Thus,  in  these  patients,  it  seems  likely  that  HFV  could  be  used 
without  undue  risk  of  barotrauma.  The  failure  of  mean  central  airway 
pressure  to  reflect  peripheral  pressures  in  the  lung  implies  that  mean 
airway  pressure  monitoring  alone  may  not  ensure  safe  use  of  HFV.  However, 
until  the  clinical  settinga  in  which  HFV  can  be  safely  employed  are 
defined  more  precisely,  lung  volume  or  another  index  of  alveolar  pressure 
should  be  monitored  to  minimize  the  risks  of  adverse  side  effects. 

E.  Military  Implication* 

These  findings  may  be  of  value  in  designing  strategies  to  deal  with 
respiratory  failure  in  field  casualties  in  the  event  of  nerve  gas 
poisoning.  Specifically,  if  the  nerve  gas  poisoning  results  in  a  loss  of 
pulmonary  eampllanoe  in  excess  of  an  Increase  in  resistance,  then  HFV  may 
be  applied  without  a  substantial  chance  of  dynamic  hyperinflation 
occurring.  In  contrast,  if  increase  in  airway  resistance  is  prominent, 
then  dynamic  hyperinflation  may  occur.  Farther,  given  knowledge  of  the 
relative  change  in  the  various  parameters  which  may  occur,  one  could 
select  possible  optimal  frequencies  and  tidal  volumes  to  use. 

III.  MODEL  STUDIES  AT  M.I.T. 

A.  Introduction 

Various  approaches  have  been  used  to  develop  methods  for  predicting 
pressure  differences  in  the  pulmonary  airways.  In  the  case  of  normal 
tidal  breathing,  theoretical  and  semi-empirical  methods  have  been 
developed  which  are  in  reasonable  agreement  with  physiologic  measurement 
(It, 15).  These  methods  of  analysis,  however,  are  typioally  confined  to 
inspiratory  flow  and  are  purposely  restricted  to  situations  analogous  to 
tidal  breathing.  In  such  cases,  the  pressure  change  is  due  to  the 
combined  effects  of  wall  shear  stress  within  the  boundary  layers  that 
develop  at  each  new  bifurcation,  and  of  fluid  acceleration  (or 
deceleration),  often  termed  *tbe  Bernoulli  effect. " 


Experimental  measurements  of  pressure  drop  across  the  airvays  are  most 
conveniently  expressed  la  dimensionless  terms  as  (15): 


l/^pur  "  fn*Re»  a*  8«onetry)  (1) 


Here  Re«ud/v  is  a  characteristic  Reynolds  number,  d  Is  the  airway 
diameter,  u  Is  the  cross-sectional  mean  velocity,  and  a«(d /2\/ui/v  Is  a 
measure  of  the  unsteadiness  of  the  flow  where  u  Is  the  osoillation 
frequency  and  v  la  the  kinematic  vlsooslty  of  the  fluid.  In  all  studies 
mentioned  thus  far,  the  dependence  upon  a  was  deemed  unimportant. 

Relatively  few  studies  have  examined  flows  that  are  intrinsically 
unsteady.  JaffTln  and  Hennessey  (16)  measured  the  pressure  drop  in  a 
model  of  the  central  airways  for  a  purely  sinusoidal  flow.  They  reported 
pressure  signals  that  possessed  both  a  nonlinear  character  and  a  phase 
shift  with  respect  to  the  flow.  Their  experiments,  however,  were  limited 
to  a  narrow  range  of  dimensionless  frequencies  (1<  cK2.66)  making  It 
difficult  to  Infer  the  general  behavior  of  unsteady  flows. 

More  recently,  Isa bey  and  Chang  (17)  reported  measurements  of  the 
total  pressure  drop  during  periodic  flow  across  a  cast  of  the  husan 
central  airways  (5  generations).  They  presented  their  results  by 
plotting  Ap/pu 2/2)  against  Re,  where  Ap  is  the  pressure  drop  at  any 
Instant  during  the  cycle  and  u  and  Re  are  the  corresponding  values  (at  the 
same  Instant  In  time)  of  tracheal  flow  velocity  and  Reynolds  number.  They 
found  that  to  completely  describe  their  data  It  was  necessary  to  Introduce 
a  new  parameter  e  «  L(du/dt)/u2.  This  parameter  can  be  rearranged  and 
written  as  (2L/d)(c^/Re)/tan  tot.  Expressed  in  this  fora,  £  is  seen  to 
be  dimensionless  time  and,  as  such,  clearly  establishes  that  the  phase  of 
oscillation  Is  of  inportanoe  In  understanding  the  physios  of  the  flow. 
However,  it  Is  possible  to  completely  describe  the  dimensionless  pressure 
drop  using  a  function  of  the  fora  of  Equation  (1)  without  the  need  for 
Introducing  e« 

In  thla  report  we  aubmlt  experimentally  derived  generalised  laws  that 
could  be  uaed  to  predlot  Ap  across  either  the  entire  lung  or  individual 
airways  over  a  wide  range  of  osolllatory  flow  conditions.  To  accomplish 
this  we  use  a  methodology  for  decomposing  the  pressure  signal  using 
fourler  analysis  and  apply  this  methodology  to  measurements  of  pressure 
drop  at  various  sites  of  a  rigid  four  generation  symmetric  model  of  the 
airways.  This  analysis  results  In  universal  curves  oapsble  of  describing 
all  our  data.  He  then  use  these  general  curves  to  predlot  the  pressure 
drop  In  a  dog  lung  model  and  compare  the  predictions  to  experimental  data. 


Th«  motivation  for  this  study  la  derived  froa  the  need  to  better 
understand  the  pressure  excursions  associated  with  high  frequency,  low 
voluae  ventilation  (HFV).  Knowledge  of  the  pressures  produced  during  HFV 
could  be  important  for  assessing  the  dangers  of  barotrauma  and  lupalred 
cardiac  perforuance  prevalent  in  other  foras  of  ventilatory  assist,  and  in 
the  design  of  high  frequency  ventilators.  These  pressure  excursions  also 
influence  the  distribution  of  Inspired  gas  and  the  aaxlaua  allowable 
voluae  flow  rate. 

B.  Methods 

1.  Pets  aamliltlon.  A11  aeasureaents  were  aade  on  a  four 
generation  (with  the  trachea  being  generation  zero)  syanetrlcally 
branching  network  asseabled  frca  eonaereial  nalgene  'I*  oonnectors  (Figure 
6).  Bach  generation  had  an  Internal  disaster  of  1  ca,  a 
length-to-dlaaeter  ratio  of  3.5  and  a  branching  angle  of  70  deques. 

The  systea  was  excited  at  the  'trachea'  with  sinusoidal  oscillations 
produced  by  a  not or  and  piston.  At  the  'alveolar'  end,  each  of  the 
sixteen  terainal  branches  was  connected  to  a  30  oa  long  pieoe  of  tygon 
tubing  which  was  open  to  atacsphere.  The  terainal  branches  were  extended 
in  this  way  to  produce  a  sore  uni fora  flow  distribution  In  the  branching 
network.  Since  these  extensions  posed  a  large  overall  lapedance  as 
cob pared  to  Individual  branching  sections  they  helped  to  unlforaly 
distribute  the  flow  to  all  brandies,  negating  the  tendency  for  flow  to  be 
preferentially  directed  down  the  aedlal  pathways  (18)  and  alnlalzlng  the 
effeet  of  slight  construction  Inaocvacies.  Thus  we  will  assuae  a 
synaetrlc  partitioning  of  flow  between  sibling  branches  in  a  given 
generation.  Evidence  to  support  this  assuaptlon  Is  obtained  by  a 
comparison  of  pressure  drops  at  analogous  aedlal  and  lateral  locations  and 
will  be  discussed  later. 

Two  parameters  were  recorded  during  each  run:  the  motion  of  the 
piston  (from  which  flow  rates  were  deduced)  and  the  pressure  difference 
between  successive  generations.  Differential  pressures  were  measured 
(using  a  Tall dyne  MP45  transducer  with  a  2  ca  HgO  diaphragm)  at  a  total 
of  six  aeaauraaent  locations  as  shown  In  Figure  6,  naaely: 


first  generation  —  between  stations  (1)  and  (1,2) 

second  generation  lateral  —  between  stations  (1,2)  and  (1,3) 

third  generation  lateral  —  between  stations  (1,3)  and  (1,4) 

third  generation  aedlal  —  between  stations  (2,3)  and  (2,4) 

fourth  generation  lateral  —  between  stations  (1,4)  and  (1,5) 

fourth  generation  aedlal  —  between  stations  (2,4)  and  (2,5) 


Pressure  taps  were  approximately  1  n  in  internal  diameter  and  were 
mounted  flat  with  the  wall,  perpendicular  to  the  plane  of  the 
bifurcation.  The  motion  of  the  piston  was  measured  with  a  displacement 
transducer.  Instantaneous  values  of  these  two  parameters  were  digitally 
processed  using  a  mini- computer.  Each  measurement  was  repeated  over  20  to 
100  cycles  and  an  ensemble  average  was  calculated  and  stored  for  later 
analysis. 

To  ascertain  the  effect  of  secondary  motions  we  also  examined  the 
pressure  differences  at  a  number  of  points  around  the  circumference  at 
station  (1,2)  in  some  experiments. 

2.  Data  reduction.  The  pressure  difference  between  any  two  measurement 
locations  (e.g. ,  (1)  and  (2)  in  Figure  7)  is  caused  by  three  effects: 
instantaneous  acceleration  of  the  mass  of  fluid  oontained  within  the 
control  volume  between  (1)  and  (2),  viscous  dissipation  due  to  friction, 
and  the  change  in  dynamic  head  as  the  fluid  passes  between  the  perent  and 
child  vessels.  This  can  be  expressed  in  terms  of  an  energy  conservation 
equation  for  the  control  volume  of  Figure  7: 


P2  -  P!  -  -py<3?>  +  s  +  fes^  -  e2q>  (2) 


where, 

A  A 

P2  and  Pi  are  weighted  average  values  of  pressure  at  stations  (1)  and 
(2)  defined  by  $  *  /pu.dA/Q 


p  is  the  density 
Q  is  the  flow  rate,  Q  ■  / udi 

D  is  the  total  rate  at  which  mechanical  energy  is  dissipated  by 
vlsooslty  in  the  region  between  (1)  and  (2) 

u  is  the  cross-sectional  average  velocity,  u  «  /udA/A 

0  is  the  ratio  of  actual  energy  flux  to  that  associated  with  the  same 
flow  rate,  with  a  uniform  velocity  profile  and  Is  defined  as:  B  *f 
u3dA/Au3.  Similarly,  Y  represents  the  ratio  of  aotual  momentum  flux 
to  that  for  a  uniform  velocity  flow,  Y  ■  /u2dA/Au2.  Both  have  values 
of  unity  for  a  flat  velocity  profile  while  for  a  parabolic  profile, 3  «2 
and  y*4/3. 


In  traditional  linear  analyses  of  snail  amplitude  oscillatory  notion, 
the  first  tarn  is  due  to  the  fluid  inertia  (I),  the  seoond  is  the  flow 
resistance  (ft)*  and  the  third  is  cnltted  since  it  represents  a  nonlinear 
effect.  For  sinusoidal  flow*  the  first  tern  on  the  right  hand  side  of 
Equation  (2),  being  proportional  to  dQ/dt*  appears  at  the  excitation  or 
"fundamental"  frequency.  The  last  tern,  being  proportional  to  the  square 
of  the  Telocity,  appears  at  tviee  the  excitation  frequency,  as  a  "first 
hamonlc. "  As  for  the  dissipation  tern  D/Q,  there  is  no  a  priori 
restriction  on  its  frequency  oontent.  If  dissipation  is  nalnly  confined 
to  a  Tiscous  boundary  layer  it  would  appear  at  the  fundanental  fTequancy, 
while  dissipation  due  to  secondary  flows  would  hare  both  a  fundanental  and 
first  barnonlo  contribution.  Therefore,  it  is  possible  to  obtain 
inforaation  concerning  the  contribution  of  each  effect  by  reducing  each 
signal  to  a  set  of  anplltudes  and  phases  corresponding  to  the  various 
hamonlea: 


Ad  ■  xj  [  p  sin(nu)t  +  e  )]  (3) 

n  n 

n»l 


where  pn  and  are  the  anplltude  and  phase  of  the  nth  harmonic  of  the 
pressure  signal. 

In  steady  flows,  the  dynanlc  head  (pu^/2)  is  conventionally  used  to 
nake  pressure  dimensionless;  this  choice  is  by  no  means  unique.  In 
particular  since  the  fundamental  represents  both  fluid  acceleration  and 
viscous  dissipation,  it  is  of  Interest  to  study  the  ratio, 


*  P2  “  Pj/pLUu 


(4) 


where  0  is  the  local  veloolty  amplitude  in  the  parent  branch,  L  *  Lj  ♦ 
1*2/2  and  Li  and  L2  are  the  lengths  of  the  parent  and  sibling 
branches,  respectively.  The  quantity  PLOW  is  the  amplitude  of  PL(dO/dt), 
the  pressure  drop  due  to  fluid  acceleration  alone  if  the  veloolty  profile 
were  blunt.  Henoe,  a  value  of  F  close  to  unity  would  indicate  the 
dcmlnanoe  of  inertial  effects  over  viscous  dissipation. 

It  is  also  of  Interest  to  compare  the  fundamental  to  the  predictions 
of  laminar  oscillatory  flow  is  a  straight  pipe  (19)  applied  separately  to 
eaoh  of  the  child  and  parent  tubes; 


Ap 

Kvom 


p2  ■  P1  "  I  PV1WL/M^0(a)]  sin(wt  -  j  -  e10(a)) 


(5) 


The  paraMters  M' iq(oO  and  e^g(a)  era  the  mm  as  tboM  ustd  by 
Woaeralay  (19)  and  are  graphically  praaaated  ia  Figure  8.  For  straight 
tubas*  M1 10  approaches  unity  and  ejo  approaches  aero  as  iaertla  coses 
to  dominate  viscous  dissipation  as  shown  ia  Figure  8  for  high  a. 

In  a  aiailar  fashion*  a  non-dlaenslonal  Manure  of  the  first  harmonic 
can  be  introduced  in  the  fora  of  the  ratios 


amplitude  of  the  "1st  harmonic"  in  the  pressure  signal 

C,  - - (6) 

amplitude  of  the  "1st  hamonic"  for  an  lnvlacld  stress 
with  blunt  eelocity  profiles. 

is  defined  above,  is  generally  referred  to  as  a  correction  factor. 

In  a  single  bifurcation  with  an  area  ratio  of  2  and  a  cross-sectional  Man 
velocity  in  the  parent  tube  of  anplltude  Oj*  the  den os ins tor  of  Equation 
(3)  would  be  equal  to  (3/t)(PD^2/2). 

C.  Results 

Experiments  were  conducted  over  frequencies  ranging  from  0.15  to  20  Hz 
and  tidal  volumes  in  the  range  of  3  to  77  ml.  This  corresponds  to  a  range 
of  Reynolds  numbers  (Re  «p0d/v  where  0  is  the  cross-MCtlonal  Man 
velocity  amplitude)  in  the  most  immediate  parent  brands  between  50  and 
30,000  and  dlmnalonleas  frequencies  (a)  ranging  from  1  to  15.  A 
representative  data  set  for  two  different  values  of  frequency  and 
"tracheal"  flow  rate  is  shown  in  Figures  9i  end  9B. 

Over  the  range  of  experlMnta,  we  observed  generatlon-to-generation 
pressure  differentials  ranging  between  0.001  and  5  cm  HgO.  it  was  not 
possible  to  directly  calibrate  our  transducer  for  pressures  as  low  m 
0.0 Oil  cm  H2O.  Therefore,  we  mde  several  Maauresents  in  a  long 
straight  tube  under  conditions  in  which  we  would  expect  VoMr slay's 
prediction  to  be  valid,  it  the  lowest  pressures,  the  Masurement  was 
nearly  509  below  the  theoretical  prediction  (Equation  5).  Further,  there 
was  a  lack  of  linearity  in  the  response  of  the  pressure  transducer  in  the 
range  of  pressures  below  0.1  om  HjO.  Since  we  were  unable  to  correct 
for  this  error  in  a  reliable  Mnner,  we  anticipate  errors  approaching  509 
at  the  lowest  pressures,  thoae  tests  having  both  small  a  and  small  Re. 
However,  this  also  la  the  range  where  fluid  acceleration  dominates  the 
pressure  signal  and  la  therefore  of  lesmr  interest. 

is  noted  earlier,  we  assumed,  for  the  purpoM  of  determining  flow 
amplitude  in  each  vessel,  that  the  flow  was  uniformly  distributed.  To  the 
extent  that  this  is  true,  we  would  expect  to  find  pressure  signals  of 
equal  magnitude  at  medial  and  lateral  positions  within  the  same 
generation.  We  show,  in  Figure  9,  exmplea  of  two  extreme  eases,  one  with 
nearly  identical  pressure  signals,  the  other  with  the  greatest  degree  of 
non-uniformity  we  observed.  There  are  significant  differences  at  the 
bitter  Reynolds  numbers  tested  in  generations  2  to  4,  but  these  tests 
represent  a  relatively  small  fraction  of  our  Mssurcments.  Furthermore, 
differences  between  the  lateral  and  Mdlal  measurements  are  prlaarlly 
associated  with  the  first  harmonic  of  the  pressure  signal.  The  maximum 
difference  between  two  first  harmonlo  signals  for  the  same  flow  condition 
was  309  while  the  average  difference  was  about  59. 
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The  signals  shown  in  Figure  9  and  noat  of  our  me  a* ur  aments  vers 
obtained  from  the  differential  between  top  pressure  taps.  To  obtain  a 
cross-section,  ve  examined  circumferential  pressure  differences  st  a 
station  in  sever al  ezperinents.  Ve  found  the  pressure  differential 
between  the  inner  and  outer  wall  of  the  nodel  to  be  greater  than  the 
pressure  differential  between  the  top  and  outer  location  but  neither 
difference  was  ever  greater  than  265  of  the  generation- to-generation 
differential.  Therefore,  ve  expect  this  figure  to  be  an  upper  satinets  of 
the  error  between  neasured  pressure  and  naan  cross-sectional  values  used 
in  Equation  (2). 

Fourier  analysis  of  the  data  revealed  that  the  pressure  signals  were 
aalnly  composed  of  a  "fundamental”  and  a  "first  harmonic,"  with  the  sum  of 
the  higher  harmonics  at  least  an  order  of  Magnitude  smaller.  While  the 
fundamental  was  typically  greater  than  the  first  harmonic,  the  difference 
diminished  as  Re  increased.  At  the  highest  Reynolds  numbers  tested,  the 
first  harmonic  was  at  tines  larger  than  the  fundamental.  Plots  of  the 
magnitude  (expressed  as  the  dimensionless  ratio  defined  in  Equation  A)  and 
phase  of  the  fundamental  are  provided  in  Figures  10  and  11,  for  the  three 
values  of  a  less  than  12.  Magnitudes  of  the  first  harmonic  as  represented 
by  the  correction  factor  Cj  (Equation  6)  are  shown  as  a  function  of  a  in 
Figure  12.  There  was  no  clear  dependence  of  on  Re  although  values 
for  Cj  in  generations  two  and  three  were  generally  greater  than  in  one 
and  four. 

D.  Discussion 

1.  general  description  of  the  results  for  alpha  leas  thin  12.  Our 
pressure  recordings  resemble  those  of  JaffTin  and  Hennessey  (16)  and 
confirm  their  observation  of  a  nonlinear  pressure-flow  relationship  and 
the  existence  of  a  phase  shift  between  pressure  and  flow.  Over  much  of 
the  range  of  our  experiments,  the  pressure  traces  show  clear  evidence  of 
the  presence  of  higher  harmonics.  Zt  must  be  recognized,  however,  that 
the  total  chlld-to- parent  area  ratio  in  our  hardware  model  was  2  while  in 
the  lung  the  changes  in  area  are  not  so  severe  (area  ratios  are  about 
1.2).  Therefore,  the  effeot  of  higher  harmonics  is  accentuated  in  our 
results. 

The  normalised  amplitude  of  the  "fundamental"  component  of  the 
pressure  signal  is  given  in  Figure  10  as  a  function  of  Reynolds  number  for 
three  values  of  o .  The  data  on  this  curve  are  not  restricted  to  any 
particular  measurement  location  or  generation  number  and  contain 
measurmsents  from  each  of  the  six  different  measurement  sites.  At  higher 
Reynolds  numbers,  some  differences  are  evident  between  the  medial  and 
lateral  measurements. 

Figure  10  shows  that  for  low  Reynolds  numbers  most  of  the  pressure 
slgnsl  is  due  to  fluid  sooeleratlon.  For  esoh  value  of  a ,  however,  the 
results  exhibit  a  different  character  beyond  a  particular  value  of 
Reynolds  number  which  turns  out  to  be  approximately  200  tines  the 
corresponding  value  of  ou  This  observation  is  of  particular  interest  in 
light  of  the  faot  that  for  oscillatory  flow  in  a  long  straight  pips,  onset 
of  turbulence  has  been  reported  at  a  Re/a  of  (200-500)  (20,21).  Rote  that 
when  a  *A,  the  break  oocurs  at  a  Reynolds  number  of  only  800. 
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For  Reynold*  nuabers  above  the  break  ^Pfund  (the  magnitude  of  the 
fundamental)  la  oonslderably  higher.  Thin  additional  pressure  drop  could 
be  due  to  an  inereaalng  nonuniformity  In  the  velocity  profile,  or  a 
dramatic  Increase  in  dissipation.  The  latter  Is  more  likely  and  oan 
probably  be  attributed  to  the  onset  of  turbulenoe  as  suggested  above. 

As  shown  In  Figure  11,  the  phase  of  the  fundamental  component  of 
pressure  relative  to  the  mean  flow  agrees  reasonably  well  with  values  from 
Vonersley*s  solution  for  a  straight  tube.  Bowever,  it  oan  also  be  seen  In 
Figure  11  that  the  phase  difference  follows  a  general  trend  even  for 
Re/a  >200  although  one  which  is  no  longer  oonslatent  with  the  laminar  flow 
solution  In  a  straight  tube. 

Based  on  the  observation  that,  at  high  a,  transition  Is  governed  by 
the  paramstsr  Re/a,  we  replotted,  In  Figure  13,  P  against  Re/a,  which  is 
the  Reynolds  number  based  on  boundary  layer  thickness.  The  plot  suggests 
strongly  that  for  4<  a<i2,  the  point  of  transition  Is  determined  primarily 
by  the  parameter  Re/a. 

Data  for  the  range  Re/ 00200  are  best  described  by  the  correlation: 


Ap  -  (0.20)Re0,62/a1,35  for  a  <  12  (7) 


This  correlation  was  obtained  from  a  double  least-squares  curve-fit  of  the 
data  with  a  good  fit  (r2  *  0.64). 

Thus,  while  It  appears  that  a  transition  oceurs  at  Re/ a  >200.  the 
value  of  P  above  this  transition  varies  approximately  as  (R*/a2)®*°. 

Note  inclusion  of  a2/Re  or  the  Strouhal  number  which  is  the  ratio  of 
temporal  acceleration  (proportional  to  3u/dt)  to  eonveotlve  acceleration 
(proportional  to  u3u/3x)  (6).  This  relationship  suggests  the  surprising 
result  that  viscous  effects  are  unimportant  in  determining  the  magnitude 
of  Ap/pu2  despite  the  anticipated  lnportanoe  of  viscous  dissipation. 

This  result  has  a  crude  analogue,  however,  in  steady  turbulent  flow  for 
which  Ap/pu2  Is  very  weakly  dependent  on  Reynolds  number,  thus 
signifying  that  viscosity  plays  only  a  small  role  in  establishing  the 
level  of  viscous  dissipation. 

2.  Results  for  a>l2.  Based  on  studies  of  oscillatory  flow  In  curved 
pipes  we  had  reason  to  expeot  more  complicated  flow  patterns  and  therefore 
a  different  soaling  pressure  for  a>i2.  For  example,  in  oscillatory  flow 
In  a  curved  pipe,  Lyne  (22)  has  developed  a  theory  which  predicts 
transition  from  a  two-oell  to  four- cell  seoondary  flow  vortex  pattern  at 
a«  12.9.  Likewise,  Mullin  and  Granted  (23)  observed  experimentally  the 
appearance  of  two  additional  vortloes  and  the  development  of  the  four 
vortex  system  at  a  ■  11.0.  Anticipating  qualitatively  similar  effects  In 
these  studies,  we  made  a  separate  plot  for  F  vs  Re/a  for  data  points  with 
a>l2  (Figure  14).  In  this  range,  the  pressure-flow  relationship  Is  more 
uniform  throughout  the  range  of  Re  with  no  abrupt  transition  between  the 
low  and  high  Re  range.  Unfortunately,  we  were  confined  to  a  rather  narrow 
range  of  a(iess  than  12),  and  therefore  a  double  least-squares  curve  fit 
of  the  data  points  In  this  range  was  not  statistically  meaningful. 


3.  am  ilfanUHa  lot— IliatifllU  For  the  ukt  Of  comparison  with 
existing  studies*  which  almost  invsrlsbly  hare  normalized  Ap  using  the 
quantity  (P  02/2),  we  hare  Included  in  Figure  15  a  plot  of  APfuad^P 
0 2/2) vs.  Be  for  a<12.  The  aolld  lines  represent  double  least-square 
curve  fits  when  the  deta  are  divided  into  two  categories:  Be/a  <100  and 
Be/a  >200.  For  Be/a  <200  the  data  exhibit  the  fora 


(8) 

where  the  Strouhal  nuaber  (Be/a2)  again  appears.  Similarly,  in  the 
Be/a>200  range  all  the  regression  lines  have  the  fora* 


*w>vi/2  -  <3><£)-4 

(9) 

Both  expressions  are  consistent  with  the  correlations  given  above. 

4.  First  Hartmle.  As  aentloned  earlier*  the  fundaaental  and  first 
hsraonlcs  of  the  pressure  drop  were  the  major  signals.  The  first  haraonle 
is  ea bodied  in  the  oorreotlon  factor  (see  Equation  7)  which  acoounts 
both  for  convective  acceleration  and  for  the  contribution  of  dissipation 
to  the  first  haraonle. 

Ve  expect  the  velocity  profiles  to  be  relatively  insensitive  to 
changes  in  Reynolds  nuaber  (as  found  in  steady  flow  through  a  brandling 
network  (24),  but  to  becoae  no re  blunted  as  the  dlaenslonless  frequency 
Increases  due  to  the  reduced  tlae  for  boundary  layer  growth.  Accordingly, 
showed  no  clear  dependence  on  Be*  but  exhibited  a  definite  tendency 
to  Increase  with  increasing  a  as  deaonstrated  in  Figure  13.  This  finding 
suggests  that  dissipation  contributes  only  to  the  first  haraonle  and  that 
the  dependence  of  on  a  apparently  Is  due  to  a  gradual  rise  in  the 
effects  of  dissipation.  This  tendency  for  C]  to  Increase*  especially 
for  a>12  nay  be  essoclated  with  the  developoent  of  a  sore  coaplex 
secondary  velocity  pattern  as  discussed  above.  A  double  least-squares 
curve  fit  yields  the  following  foras  for  C^: 


1.41  for  a  <  10 
-1.93  +  0.37  a  for  a  >  10 

(10) 

Our  findings  that  the  pressure  signal  can  be  adequately  described  by 
the  first  two  baraonlcs  suggests  that  the  set bod  of  analysis  proposed 
here*  while  soaeWhat  comparable  to  that  of  Chang  and  Isabey  (17)*  baa  soaa 
practical  advantages.  The  differences  between  the  two  approaches  can  be 
illustrated  by  attempting  to  predict  the  presaure  signal  (as  aeasured  in 
the  trachea,  for  example)  associated  with  a  particular  flow  condition. 


Using  tbs  method  of  four lsr  summation,  it  oan  bs  obovn  that  it  is 
suffldsnt  to  consider  tbs  first  two  terms  of  tbs  ssriss  in  Equation  (3) 
sines  all  blgbsr  harmonics  ars  small  by  oompsrlson.  Both  tbs  coefficients 
of  this  ssriss  and  tbs  phasss  arc  simple  functions  of  tbs  tun 
dlmensionlsss  groups,  la  and  a  (which  can  bs  calculated  from  tbs  airway 
glows  try  sad  local  flow  conditions)  and  ars  given  by  tbs  correlation  for 
pressure  amplitude  (Equation  7)  and  Cj  (Equation  10)  and  by  assuming 
that  the  fundamental  is  of  the  same  phase  as  predicted  for  a  straight  tube 
and  the  first  harmonic  is  in  phase  with  the  square  of  the  time-Yarying 
flow  rate.  These  coefficients  can  be  used  in  Equation  3  to  yield  tbs 
time-varying  pressure  amplitude  and  phase  across  each  branch  which  can 
tben  be  summed  over  all  the  individual  airways  to  obtain  the  overall 
pressure  difference.  With  Xsabey  and  Chang's  approach  (17)»  it  is 
necessary  to  perform  a  sufficient  number  of  experiments  so  that  a  general 
expression  for  dimensionless  pressure  P  (or  Ap/pu2)  oould  be  obtained. 

The  dependence  of  P  on  e  will  describe  the  time-variation  of  the  pressure 
during  s  single  cycle. 

5.  Caanarlson  to  Physiologic  Experiments  on  Doss.  As  noted  above,  these 
results  are  subject  to  error  due  to  differences  between  the  model  and  the 
pulmonary  airways  as  they  exist  in  vivo.  First,  the  area  ratio  between 
parent  and  child  veasels  was  higher  in  our  experiments  than  in  the  lung. 
This  would  tend  to  overemphasise  the  effects  of  oonveetlve  acceleration 
and,  la  combination  with  the  sharp  corners  within  the  bifurcation,  would 
promote  separation  and  the  onset  of  turbulence.  Second,  in  situations  of 
high  flow  rate,  compliance  may  play  a  significant  role  by  producing 
changes  in  airway  geometry. 


To  test  the  applicability  of  the  present  results  to  the  slrways,  and 
in  particular  to  examine  the  improvement  obtained  by  this  approach  over 
the  linear  theory,  a  simple  model  was  developed  to  estimate  total  airway 
resistance  in  a  dog  lung  under  typloal  EFT  conditions.  In  this  section  we 
compare  the  predictions  of  this  model  to  the  measurements  of  total 
tranapulmonary  pressure  drop  in  five  dogs  neasired  st  the  Harvard  School 
of  Publlo  Health. 


The  theoretleel  cel eolations  were  based  on  the  10  generation  symmetric 
model  of  the  dog's  bronohlsl  tree  of  Horsfleld  and  Cunsing  (25).  Ve 
assumed  the  airways  to  be  rigid  and  the  gas  to  be  incompressible.  The 
effects  of  lnertls  and  viscous  losses  due  to  the  motion  of  the  gas  ware 
modeled  according  to  the  expressions  given  above. 

When  a  specific  tidal  volume  and  frequency  are  imposed  at  the  trachea, 
the  assumption  of  symmetric  portioning  of  the  flow  between  the  sibling 
branches  in  s  given  generation  allows  one  to  calculate  local  values  of 
Reynolds  number.  Once  Re  and  a  in  each  generation  are  known,  the 
correlations  obtained  in  this  study  can  be  used  to  determine  the 
magnitudes  of  the  fundamental  and  first  harmonic  of  the  pressure 
difference  across  that  generation.  The  fundamental  was  assumed  to  have 
the  same  phase  as  would  exist  in  a  straight  tube  at  the  corresponding 
value  of  a,  while  the  first  harmonics  were  assumed  in  phase  with  the 
flow.  Calculation  of  total  (tracheal-to-alveolar)  pressure  amplitude  was 
performed  by  adding  up  the  contributions  from  the  fundamentsl  end  first 
harmonic  oomponents  at  each  of  the  ten  generations  of  the  model.  In  all 
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ouu  considered,  the  onrill  contribution  of  tbo  first  barnonlc,  ones 
msm4  our  all  tbs  too  generations  of  the  not  work,  was  negligible.  lbs 
cannula  tl  vs  prsaaurs  amplitude  duo  to  both  components  was  negligible 
beyond  this  point. 

Vs  also  incorporated  tbs  sffsot  of  lone  parenchyma  and  obsst  wall 
lapedanoe.  These  sffsota  wars  ooablnsd  to  fora  a  tarn  Inal  lnpsdanos, 
zsad  *  *sff  ♦  wCeff*  tbs  lnsrtanos  asaoolatsd  with  tbs  chsst 
wall  was  assuaad  to  ba  snail.  Taluss  of  H*ff  and  C*ff  wars  sstiaatsd 
by  comparing  total  prsaaurs  drop  (dus  to  air  notions  only)  to  nsasursd 
prsssurss  at  tbs  two  lowsat  valuss  of  fFequeney  (2  and  t  Ha)  and 
attributing  tbs  sntlrs  dlffsrsnos  to  this  calculation  gave  a 

valus  for  tbs  effective  oonplianos  of  lung  tissue  and  thorax  C^ff  « 
0.023/L/onB20,  and  an  sffsetlvs  rsslstanos  of  lung  tissue  and  thorax 
B«ff  *  2.5  on  BjO/l/sso,  valuss  which  are  eonparabls  to  the 
nsasursnsnts  of  Crossfill  and  Viddieoabs  (26). 

Using  these  valuss  of  R*ff  and  C^ff  the  total  transpulaonary 
pressure  drop  in  the  dog  lung  was  cooputed  at  each  value  of  frsqusnoy  and 
tidal  voluae  used  in  the  sxpsrinsnts.  lbs  predictions  of  tbs  theoretical 
nodsl  and  its  comparison  to  in  vivo  nsasursnsnts  are  shown  in  Figure  16. 
Further,  In  order  to  highlight  the  findings  of  this  study,  the 
contribution  to  the  total  pressure  anplitude  due  to  parenohynal  and  obest 
wall  effects  alone  is  also  shown.  The  graph  shows  renarkable  agreensnt 
between  the  present  nodsl  and  the  sxpsrinsnts  and  also  dsnonstrates  that 
while  tbs  linear  predictions  (treating  the  airways  as  individual  straight 
tube  segMnta)  uy  be  adequate  for  very  low  tidal  voluae  frequency 
ooablnatlons  it  underestimates  the  pressure  drop  over  such  of  the 
practical  range  of  HFV.  The  nodsl  has  obvious  utility  in  calculating 
pressure  drops  that  nay  occur  during  HFV  in  the  airways. 


IV.  Floras  LEGENDS 


Figure  1:  Clroult  for  high  frequency  ventilation  (HFV). 

Figure  2:  Lung  volume  Iootum  during  high  froquoncy  ventilation  as  a 
function  of  froquoncy  at  various  tidal  volunes.  Each  syabol  ropresonts  a 
tidal  voluae:  solid  circle,  50  al;  solid  square,  100  al;  and  open  circle 
150  al. 

Figure  3:  Lung  voluae  increase  plotted  ss  a  function  of  tracheal  flow 
aaplltude  (f  x  V  T*)  during  high  frequency  ventilation.  Each  ayabol 
represents  a  tidal  voluae:  solid  circle,  50  al;  solid  square,  100  al;  and 
open  circle,  150  al. 

Figure  V:  Lung  voluae  Increase  during  high  frequency  ventilation  (HFV) 
plotted  as  a  function  of  respiratory  aystea  relaxation  pressure  (F  stat*) 
after  HFV  in  aix  patients.  Relaxation  preaaures  were  not  obtained  in 
patient  8  and  no  change  in  lung  voluae  was  noted  in  patient  1. 

Figure  5:  Panel  A  displays  a  recording  of  airway  preaaure  during  high 
frequency  ventilation  (HFV).  The  phasic  tracing  was  apparent  only  when 
the  electrloal  averaging  olroult  waa  switched  off.  Panel  B  shows  lung 
voluae  during  conventional  ventilation  and  HFV.  At  "a"  the  tracheal  tube 
is  oonnected  to  the  weighted  splroaeter,  and  lung  voluae  increases.  At 
•b"  HFV  begins  and  a  further  increase  in  lung  voluae  occurs. 

Figure  6.  Branching  network  used  in  the  experiments.  Heasureaent 
locations  are  indicated  by  nuabers  la  parentheaes.  Where  two  nuabers  are 
given,  the  first  indloates  either  aedial  (2)  or  lateral  (1)  location;  the 
seoond  indloates  generation  nuaber  with  the  "trachea"  denoted  as 
generation  1.  Long  tubes  added  at  the  "alveolar"  end  are  not  shown. 

Figure  7.  Control  voluae  used  in  the  derivation  of  the  energy 
conservation  equation  (5)  which  i»  also  given.  The  different  eoaponents 
which  ccaprlse  Ap  are  shown  to  contribute  either  to  the  fundanental  or 
first  haraonlc  la  the  fourler  series  representation. 

Figure  8.  The  paraaeters  ejg  and  Mf ig  froa  the  solution  for  Ap  given 
by  Woaersley  (1955).  Mote  that  over  aucb  of  the  range  of  these 
experiments  (*<  <1A)  H'lg  is  olose  to  unity  and  e^g  i»  less  than  3°'. 

M'io  represents  the  ratio  pGLu/Apygg. 

Figure  9*  Coaparlsons  between  the  aeasured  pressure  signal  between  a 
aedial  (solid  line)  and  lateral  (dashed  line)  location  within  the  saae 
generation  for  the  saae  flow  conditions:  (a)  Re* 136,  0*7.27  in  the  third 
generation,  (b)  Re* 1296,  0*7,23  in  the  third  generation,  (a)  and  (b) 
represent  the  olosest  and  poorest  agreeaent,  respectively,  between  the 
aedial  and  lateral  locations. 

Figure  10.  Dimensionless  pressure  plotted  against  the  Reynolds  nuaber 
for  a<12.  •,  a*A,  ■,  a*7,  ▲,  a*10.  Open  synbols  lndloete  lateral 

aeasureaents,  olosed  symbols  indicate  aedial  aeasureaents  or  aeasureaents 
In  the  first  or  seoond  generations. 


Fleur*  11.  Th*  neaaured  phase  dlff*r*oo*  between  the  fundsaental 
component  of  the  pressure  signal  end  the  flow  Telocity.  In  the  Halt  of 
high  frequency i  pressure  lends  th*  Telocity  by  90  .  Also  shown  Is  the 
phase  difference  predieted  for  oscillatory  flow  of  the  ease  diaenslonlesa 
frequency  la  a  straight,  unifora  disaster  tube  (19). 

Figure  12.  The  correction  faetcr  plotted  against  diaenslonlesa 
frequency.  Bara  represent  on*  standard  derlatlon  froa  th*  seen  at  each 
Talus  of  a. 

Figure  13.  Diaenslonlesa  pressure  saplltude  of  the  fundsaental  ooaponent 
plotted  ablest  Xe/B  for  <K12.  • ,  M,  I,  °»7»  a,  QtclO.  Open 

syabols  Indicate  lateral  aaaaursnents,  closed  syabola  Indicate  aedlal 
asasureasnts  or  aeaaursaents  la  the  first  or  second  fenerations. 

Figure  It.  Diaenslonlesa  pressure  saplltude  of  the  fundsaental  ooaponent 
plotted  against  Re/a  for  a>12.  B ,  a«l2,  •  ,  <WU.  Open  syabols 

Indicate  lateral  aeaaursaents ,  closed  syabols  Indicate  aedlal  aeasursaents 
or  asasursasnts  la  th*  first  or  seoond  generations. 

Figure  15.  Fundsaental  pressure  saplltude  Bade  dlasnslonless  by  PU2^ 
plotted  against  Reynolds  nuaber.  ▲  ,  cuio,  B ,  a*j,  • ,  oat.  Lines 
represent  best  least  squares  fit. 

Figure  16.  Ccaparlson  between  neasured  pressure  saplltude  in  th*  trachea 
of  anesthetlsed  dogs  Indicated  as  a  a* an  and  standard  derlatlon  of  each 
set  of  aeasureaenta,  and  pressures  predicted  by  these  hardware  ezperlasnts 
(solid  line),  by  aassalng  that  Vcaarsley**  solution  la  appropriate  (dashed 
line),  and  by  neglecting  airway  resistance  entirely  (dot-dash  line). 
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Table  1 

Characteristics  of  Patients 


Patient  Sex  Age 

(yr) 

1  F  64 

2  F  55 

3  F  30 

4  F  34 

5  M  50 

6  F  56 

7  F  64 

75 


41.6 

Bronchiectasis, 
Pulmonary  Fibrosis 

40.1 

Amyotrophic 

Lateral 

Sclerosis 

44.9 

Cerebrovascular 

Accident 

46.1 

Cerebrovascular 

tecident 

61.2 

Anoxic 

Encephalopathy 

43.3 

Asthma,  Hepatic 
Encephalopathy 

49.9 

Syringcnelia 

49.8 

Asthma 

Static  Compliance 
(1/an  HjO) 

.012 

7  .040 

6  .045 

6  .059 

8  .037 

8  .032 

4  .063 

7  .021 


Tracheal 

Weight  Diagnosis  Tube 

(kg)  (ran  i.d.) 

10 
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